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Many cytoskelctal biopolymers are "active," consuming energy in large quantities. In this Letter, 
we identify a fundamental difference between active polymers and passive, equilibrium polymers: 
for equal mean lengths, active polymers can reorganize faster than equilibrium polymers. We show 
that equilibrium polymers are intrinsically limited to linear scaling between mean lifetime and mean 
length, MFPT ~ (L), by analogy to 1-d Potts models. By contrast, we present a simple active- 
polymer model that improves upon this scaling, such that MFPT ~ (L)^^^. Since to be biologically 
useful, structural biopolymers must typically be many monomers long, yet respond dynamically to 
the needs of the cell, the difference in reorganization kinetics may help to justify active polymers' 
greater energy cost. PACS numbers: 87.10.Ed, 87.16.ad, 87.16.Ln 



Cytoskeletal polymers play a key role in cellular repro- 
duction, locomotion, and transport [1-3]. Biopolymers 
like actin filaments and microtubules in eukaryotes and 
FtsZ, MreB, and ParM in prokaryotes grow by accumu- 
lating monomers bound to the nucleotide triphosphates 
ATP or GTP. The monomers hydrolyze these triphos- 
phates to the diphosphates ADP or GDP, consuming en- 
ergy in an irreversible process and inducing conforma- 
tional changes that destabilize the polymers. In some 
cells, cytoskeletal ATP consumption can approach 50% 
of total cellular ATP consumption [4, 5]. What advan- 
tage do active polymers offer over passive, equilibrium 
polymers to justify this costly energy expenditure? 

We highlight a fundamental difference between active 
and equilibrium polymers: active polymers can reorga- 
nize faster than equilibrium polymers. Moreover, this 
difference in reorganization times widens as mean poly- 
mer length grows. Since biological structures like mitotic 
spindles or pseudopods must reach a certain size to ac- 
complish their function, yet be quickly deconstructed and 
reorganized, this intrinsic difference may at least partly 
justify active polymers' greater energy cost. 

A large class of equilibrium models describes a poly- 
mer as an ordered sequence of monomers, each of one 
of q types (including different conformational states of 
the same molecule). Monomers can attach, detach, and 
potentially interconvert among the q types. Interactions 
between neighboring monomers {?',?' + 1} contribute free 
energy J^^^i_^.ly to the total free energy of the polymer. 
Such models can also describe polymers consisting of 
bimdles of k protofilaments, by increasing the number of 
"monomer" types to g*'. These models are generalizations 
of 1-d, g''^-state Potts models [6]. The free energy of an 
equilibrium polymer in these models scales as the poly- 
mer length L for large L, specifically ~ LAmax) where 
Amax is the largest eigenvalue of the transfer matrix [7]. 
Hence, the equilibrium distribution of polymer lengths 
will be exponential, p{L) ~ e~^^^^\ with a character- 



istic mean length (L) = -J^. Because for large L the 
free energy effectively depends on only the largest eigen- 
value Ainax, the dynamics are essentially one-dimensional 
even for polymer bundles. This means that the effective 
force — ^ ~ — Aniax is constant, generating a constant 
negative-velocity drift in the polymer length, with drift 
velocity Vd. cx — Amax- (Polymers maintain a finite equilib- 
rium distribution because this negative drift is balanced 
by diffusion and nucleation of new polymers.) Impor- 
tantly, since the polymer length drifts towards zero at 
constant negative drift velocity, starting from the nucle- 
ation length Lnucb the mean polymer lifetime scales as 
"'l"^!' cx (^)- Thus the mean polymer lifetime 

or "mean first-passage time" (MFPT) scales as the mean 
length (L) for equilibrium polymers. This linear scaling 
is a fundamental limit for an equilibrium polymer. In or- 
der to improve upon it, a biological system must employ 
active or out-of-equilibrium processes. As an example, 
we present a simple active-polymer model based on mi- 
crotubule dynamics that yields MFPT (L)^^"^- 

Microtubule growth and disassembly dynamics have 
been well-studied [8-10]. In microtubules (and ParM), 
GTP hydrolysis leads to stochastic rapid disassembly of 
the entire polymer in a process called dynamic instabil- 
ity; the classic experimental results ]11] are reviewed in 
[1]. Recent detailed models aim to explain specific as- 
pects of the experimental data [12-16]. We consider in- 
stead a minimal microtubule model [17] that incorpo- 
rates dynamic instability. Specifically, we model an ac- 
tive polymer as an ordered sequence of monomers, each of 
which is bound either to GTP or GDP (Fig. 1). We call 
the group of GTP-bound monomers at the front of the 
polymer the "cap" and denote its size by x. We denote 
the total number of monomers (the polymer length) by 
L. GTP-bound monomers bind and unbind at the front 
end of the polymer with rates /c+ and respectively 
(Fig. lA). GTP-bound monomers at the back of the 
polymer cap undergo hydrolysis to become GDP-bound 
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Figure 1: Schematic of growth and disassembly of a model 
active polymer. (A) GTP-bound monomers (green) bind and 
unbind at the front end of the polymer. (B) The last GTP- 
bound monomer at the back of the "cap" undergoes hydrolysis 
to become an GDP-bound monomer (blue), decreasing the 
cap size by one. (C) If the cap size shrinks to zero, the polymer 
completely disassembles. 



monomers with rate fc/j (Fig. IB). If the cap size shrinks 
to zero, the polymer completely disassembles (Fig. IC). 
New polymers of length and cap size 2 are nucleated with 
rate fcnuci- We call the concentration of free GTP-bound 
monomers c and analytically treat the mean-field regime 
where c is constant, a good approximation for eukaryotic 
cells where the number of monomers is typically large 
(~ 10^). For comparison, we also consider an equilibrium 
polymer that obeys the same rules but without hydrol- 
ysis so that its length and cap size are equal. We show 
explicitly that this particular equilibrium model satisfies 
the general equilibrium scaling relation MFPT ^ (L) . 

The exact master equation for the concentration C^^x 
of polymers of length L and cap size x is 



(3) 

where the initial condition is -L, 0) = 5{x — 2)S{L — 2), 
i. e. polymers nucleate with length and cap size 2, and the 
boundary condition is p(0, -L,t) = 0, i.e. polymers with 
cap size zero disassemble. (Changing the nucleation size 
does not affect any essential results.) 

The distribution of polymer lifetimes or "first-passage 
times" (FPTs) is 
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and the mean first-passage time (MFPT) has the simple 
form MFPT = —2/ a. The steady-state polymer length 
distribution is 



-Pactivo {L) ^ \j p{x,L,t)dxdt , 



yielding the active-polymer average length 
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so that for long polymers ((i) S> 1, \a\ <C 1), one finds 
(L) ~ gD/a^, and therefore MFPT = ~2/a - (L) 
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Coarse-graining this equation leads to a continuum 
Fokker-Planck (FP) description of the probability p — 
p{x,L,t) that a single active polymer will have length L 
and cap size x at a time t after its birth: 



leading order. This sublinear scaling requires some non- 
equilibrium process, here the irreversible hydrolysis of 
monomers. Note that the details of the non-equilibrium 
model do matter to the degree of sublinearity; for ex- 
ample, a more realistic microtubule model [17] yields 
MFPT - (L)° '^ [18]. 

By comparison, in the equilibrium limit of this model 
the cap is the entire polymer. Dropping the d/dL terms 
in (2) yields an equilibrium solution that looks like (3) 
without the factor S(L — gt — 2). The equilibrium length 
distribution can then be obtained by integrating over 
polymer age: 
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where D ^ D^^ - ^ 
^{k^c + k^), and D^l — k^c + fc_ are diffusion co- 
efficients, a — fc+c — fc_ — kh is the cap drift velocity 
(a < 0), and g = fc+c — fc_ is the length drift velocity. 
The FP equation (2) describes the time-evolution of an 
individual polymer born at time (the time coordinate 
t represents polymer age). In the mean-field regime each 
polymer evolves independently once nucleated, and hence 
the fcnuci term does not appear (for details see [18]). As- 
suming g ^ —a, the effect of cap diffusion dominates the 
effect of length diffusion, and so we may neglect the Dll 



and the equilibrium-polymer average length is {L)^^^^^ ~ 
—D/a. Hence, we recover the linear scaling MFPT ~ (L) 
expected for equilibrium polymers. 

To validate these scaling relations, the master equation 
(1) was simulated using the Gillespie algorithm ]19]. The 
monomer addition rate was held fixed throughout the 
simulations. For equilibrium polymers we set fc/j = 0, 
while for active polymers for simplicity we set fc_ = 0. 
The nucleation rate fcnuci was tuned to hold the steady- 
state fraction of polymerized material constant at 75%, 
with fcnuci from 10^^ to 10~^. This is in line with experi- 
mentally measured values for microtubules ]20]. (Chang- 
ing the fraction of polymerized material to 95% does not 
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affect the qualitative results [18].) This leaves a sin- 
gle free parameter, fc_ for equilibrium polymers and kh 
for active polymers, to control the polymer length and 
MFPT. With only a single free parameter, each system 
is fully constrained by holding either the MFPT or the 
average length fixed, thus yielding a fair comparison be- 
tween the equilibrium and active polymers. 

Figure 2A shows the MFPT of the equilibrium and 
active polymers as functions of their average length. 
The data points are from simulations using 400,000 
monomers: the curves are from MFPT = — 2/a combined 
with = -D/a and (6) for (L)^,,.,,. The MFPT 

1/2 

scales ^ {L) for equilibrium polymers and ^ {L) ' for 
active polymers as expected. Hence, for the same aver- 
age length, the active polymers have much shorter mean 
lifetimes than the equilibrium polymers, and this differ- 
ence widens as average length grows. Figure 2B compares 
length distributions for the equilibrium and active poly- 
mers with the same MFPT (~10). Theoretical results 
from (5) and (7) are shown in black. Agreement between 
simulation and theory is excellent, validating our use of 
the FP equation. (For a comparison of length distribu- 
tions with the same (L), see [18].) 

What might be the biological consequences of the dif- 
ferent equilibrium and active polymer scaling relations? 
To address this question, we examine the time scale for 
large-scale spatial reorganization of structures, i.e. the 
time needed for a system to disassemble polymers at one 
site, move the material to another site, and reassemble 
new polymers. Cells often accomplish large-scale poly- 
mer reorganization in vivo by spatially regulating nu- 
cleation [21-23]. To model such regulation simply, we 
consider two spatial sites. We start simulations with nu- 
cleation occurring only at site 1 , allow the system to come 
to steady state, then switch off nucleation at site 1 and 
switch on nucleation at site 2. Monomers are assumed to 
transition between the two sites with a "diffusion" rate 
ku, while polymers do not diffuse. We define the "re- 
organization time" as the time needed after the nucle- 
ation switch for 50% of the final steady-state amount of 
polymerized material to assemble at site 2. Initially we 
assume diffusion to be fast (ku = cx)) so that a single ef- 
fective pool of free monomers is shared between the two 
sites, and then we consider the more biologically-relevant 
finite-diffusion regime. 

Figure 2C shows the reorganization time for our simple 
equilibrium and active polymers as functions of their av- 
erage length, with fast diffusion. The reorganization time 
scales ~ (L) ^ for the equilibrium polymers and ^ (L) for 
the active polymers. Hence, for a given average length, 
the active polymers reorganize faster than the equilib- 
rium polymers, and like for the MFPT the difference 
widens with increasing average length. The scaling ex- 
ponents differ from those for the MFPT because the re- 
organization time is dominated by a few very long-lived 
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Figure 2: Lifetimes versus length for equilibrium (blue) and 
active (red) polymers. (A) Mean first-passage time (MFPT) 
to polymer disassembly versus average length. The monomer 
addition rate k+ is held fixed, and time is always measured 
in units of (fc+c)~^, the mean time for monomer binding. 
Solid lines are fits to theory with slopes 1 and 0.5 for equilib- 
rium and active polymers, respectively. Some error bars are 
smaller than data-point symbols. (B) Length distributions 
with MFPT~10, with theoretical fits (black) from (5) and 
(7). The equilibrium polymers have average length 30, while 
the active polymers have average length 600. (C) Reorganiza- 
tion time (defined in text) versus average length. Solid lines 
are theoretical fits from (8) with slopes 2 and 1 for equilib- 
rium and active polymers, respectively. (D) Reorganization 
time versus monomer diffusion constant. Dashed lines show 
infinite-diffusion limits. Insets: time series showing average 
polymer number (left) and average polymer length (right) at 
site 1 versus time after nucleation is switched off. 



long polymers, whereas the MFPT is dominated by many 
short-lived short polymers. 

To understand the scaling relations in Fig. 2C, we 
estimate the reorganization time analytically as half the 
"average material age," which we define as the average age 
of the polymers in a steady-state snapshot of the system, 
weighted by the length of each polymer: 



1 age{L)P{L)LdL 

2 j^P{L)LdL 



(8) 



This average material age captures the amount of time 
an average monomer spends in an active polymer before 
it turns over, and well-approximates the reorganization 
time. Long polymers have more material than short poly- 
mers, and so the average material age is weighted by the 
length of each polymer. For our equilibrium polymers, 
age(L) = -i/a, which gives t^^,"" ~ 2D/a'^ - {if. 
Indeed, since the drift velocity a ~ 1/ (L) for any equi- 
librium polymer model as discussed above, age(L) ~ L^, 
and thus the equilibrium reorganization time scales gen- 
erally ^ {L)'^ ■ As for the MFPT, this scaling is a 
fundamental property of equilibrium polymers. To im- 
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prove upon it requires active energy-dissipation or some 
non-equilibrium process. For example, our simple ac- 
tive polymers have age(L) = L/g and hence r^eOT™ ~ 

'^^'^ 3a^D-a)^ ^ ' ^° large (L) with \a\ <C! 1, 

Next we consider the effects of slow monomer dif- 
fusion. Figure 2D shows the reorganization time ver- 
sus the monomer transition rate between sites for 
equilibrium and active polymers with the same average 
length (~90) and fraction of polymerized material (75%). 
The fast-diffusion limits, which arc reached for fcu ~ 1 
(in units of k+c), are shown with dashed lines. Strik- 
ingly, the active polymers still reorganize much faster 
than the equilibrium polymers even when slow diffusion 
limits the reorganization time. To understand this ef- 
fect, consider the equilibrium polymer dynamics after 
nucleation is switched off at site 1: polymers there dis- 
assemble stochastically and, since diffusion is slow, the 
released monomers typically rejoin other polymers at 
the same site. Hence the number of polymers at site 
1 drops rapidly while the average polymer length grows 
(Fig. 2D, insets). Eventually site 1 has only a few very 
long polymers. The equilibrium polymers remain in this 
state for a very long time, exchanging monomers with 
the free monomer pool at site 1 while diffusion slowly 
drains the pool, until the polymers finally disassemble. 
Therefore, the time needed for slow diffusion to move half 
the monomers from site 1 to site 2 is a good rough ap- 
proximation for the equilibrium reorganization time [18]. 
In contrast, active polymers do not release monomers to 
the free monomer pool except by disassembly. After the 
switch in nucleation, the number of active polymers at 
site 1 drops while the average length grows, similar to the 
equilibrium polymers. However, the few long-lived active 
polymers at site 1 then quickly accumulate and hydrolyze 
all the free monomers there, and then all the polymers 
disassemble. The time to hydrolyze all the monomers at 
site 1 . plus the time for half of those monomers to diffuse 
to site 2, is therefore a good rough approximation for the 
active reorganization time [18]. 

In summary, we find a fundamental difference between 
active and equilibrium polymers: active polymers can 
reach a fixed mean length with faster reorganization ki- 
netics than equilibrium polymers. Very generally, we 
show that equilibrium polymer lifetimes scale linearly 
with mean length. In contrast, active polymer lifetimes 
can scale sublinearly, for example as {L)^^^ in a simple 
model motivated by microtubules or as (i)*^'^ in a more 
realistic model [17, 18]. Furthermore, in our example 
the kinetic advantage of active polymers persists even 
for slow monomer diffusion. In a dynamic cellular envi- 
ronment, this kinetic advantage may help justify active 
polymers' greater energy cost. 

Our comparison of active and equilibrium polymers 
predicts one might find equilibrium polymers in biologi- 



cal contexts where polymer turnover is slow or structures 
rarely need to be reorganized. This may be the case for 
eukaryotic intermediate filaments [24] or for the bacte- 
rial homolog crescentin [25[. In addition, the existence of 
proteins like formins and profilins that accelerate actin 
polymerization suggests that kinetic regulation of active 
polymers is important to cells. Finally, although the spe- 
cific active model we consider is most closely based on 
polymers like microtubules or ParM that exhibit dynamic 
instability, our conclusions regarding accelerated kinetics 
could also relate to actin networks for which branching 
plays a role analogous to nucleation [26]. 

Our nioclel neglects many complexities of real biopoly- 
mers; clearly active polymers accomplish more than sim- 
ply reaching a certain length with a certain lifetime. We 
only suggest that fast reorganization kinetics might be 
a general (and generally desirable) feature of the active 
polymer systems that are ubiquitous in biology. 
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Here we solve the full Fokker-Planck equation defined in the main text, and justify the simplifications employed 
there. Next, we estimate the reorganization time for polymers to disassemble at one spatial site and reassemble at 
another spatial site, for slow monomer diffTision. We then compare equilibrlTini and active-polymer length distributions 
with the same average length, and also consider the result of varying the steady-state fraction of polymerized material. 
Finally, we obtain scaling relations for a more realistic model of microtubule dynamics [1], and demonstrate that an 
equilibrium two-species model satisfies our predicted equilibrium scaling relation. 

FOKKER-PLANCK DERIVATION 

Here we solve the full Fokker-Planck equation (2) defined in the main text, and justify the simplifications leading 
to (3). The Fokker-Planck equation for the probability p = p{x, L, t) of finding a polymer with length L and cap size 
x at a time t after the polymer's birth is 



dp 
di 
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LL 
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dp dp 



' dx dL dx 



'dL 



(SI) 



where D^x = \ {k+c +k- + kh), Dll = \ {k+c + and D^l = k+c+k- are diffusion coefficients, a = k-^.c—k-—kh 
is the cap drift velocity (a < 0), and g = k+c — k- is the length drift velocity. This equation does not contain a term 
proportional to fcnuci- The master equation (1) describes the time-evolution of an entire system of polymers (the time 
coordinate t represents wall-clock time), whereas the Fokker-Planck equation (2) describes the time-evolution of an 
individual polymer born at time (the time coordinate t represents polymer age). In the mean-field regime we treat, 
the growth of an individual polymer born at time is completely independent of the other polymers in the system 
(i.e. once nucleated, each polymer evolves independently). Therefore, the nucleation rate fcnuci does not explicitly 
affect the normalized distributions of polymer lengths, caps, or lifetimes. However, the nucleation rate does affect the 
steady-state amount of polymerized material and hence the free monomer concentration c, in cases where c is allowed 
to vary. 

The general solution to the Fokker-Planck equation (2) is 



p{x,L,t) = 



■- exp 



D^^{L-gt-2)'^+Di^L{x-at-2)'-^-D^i,{L-gt-2){x-at-2) 
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which satisfies both the initial condition p(x, L, 0) = 5{x — 2)S{L — 2) and the boundary condition p(0, L, t) = 0. In 
the limit Dll = Dxl = 0, this reduces exactly to (3). We note that we have solved (SI) subject to no boundary 
conditions on L. This is a valid assumption for our active polymers, which nucleate with length inuci = 2 and then 
grow at average rate g > 0, and therefore the effect of the physical boundary condition at L = is negligible. More 
formally, this assumption is valid provided g ^ —a, which implies the hydrolysis rate kh ^ —a. The probability that 
a polymer is still alive, i.e. has not disassembled, by age t is 



-Paiive(i) = / / p{x,L,t)dxdL , 
Jo Jo 

and the distribution of polymer lifetimes or "first-passage times" is 



PpPT(i) = -|Wi) = ^^ 
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Figure SI: Comparison of active-polymer length distributions for (L) ~ 132. Length distributions were calculated from 
Fokker-Planck theory using the full joint distribution of polymer length and polymer age (S6) in red and the delta-function 
approximation p{L,t) — PaiivG(t) ■ S{L — gt — 2) used in the main text (5) in green. 



which is identical to (4) in the main text. Finally, in order to calculate the steady-state polymer length distribution 
(5), we require the joint distribution p{L,t) between active polymer length L and age t. This may be obtained by 
integrating (S2) over cap sizes: 



and then the length distribution is 



p{L,t) ~ / p{x,L,t)dx, (S5) 







active I, 



(L) = / p{L,t)dt. (S6) 







In the limit D^l = D^l = employed in the main text, the joint distribution p{L,t) = Paiivc(i) • 6{L ~ gt ~ 2). 
Figure SI shows a direct comparison of the length distributions using both the full joint distribution (S6) in red and 
the delta-function approximation p{L,t) = Paiive(i) • 5{L — gt — 2) used in the main text (5) in green. The length 
distributions are nearly identical, providing a practical validation for the intuitive simplifications used to obtain (3)- (6) 
in the main text. Inserting the delta-function approximation into (S6) gives 

-Pactivc(i) = ^-Palivc ( ^ ) (S7) 
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and the mean polymer length is calculated by integration by parts: 



r ^active(L) dL J^^ p^^^^^ U)dL ^J-L PppT U)dL 



More formally, we may argue that the effect of cap diffusion dominates the effect of length diffusion. Motivated by 
the scaling relations found in the main text (and the fact that the average cap size (x) = — — ^ (L)^^^), we make 

the ansatz that we may rescale L to L — L/ (L) and x to x = x/ (L) ' . Then we obtain the non-dimensionalized 
Fokker-Planck equation: 

dp _ cfj)_ Dll d^p DxL d'^p a_ dp g_dp^ . . 

dt ~ (L) 9P + (L)2 dZ^ {Lf^ dxdl {Lf^ dx (L) QL' ^ ' 

For long polymers ((L) 3> 1) the -Dll ^^nd D^i^ terms will be much smaller than the remaining three terms, and we 
may neglect the Dll and D^l terms. (More generally, this approximation applies as long as kh,g ^ —a.) This leads 
immediately to the solution (3) in the main text. 
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POLYMER REORGANIZATION TIME 

Here we estimate the reorganization time for polymers to disassemble at one spatial site (site 1) and reassemble at 
another spatial site (site 2) following a switch in nucleation. As defined in the main text, the reorganization time is 
the time needed after the switch for half of the final steady-state amount of polymerized material to assemble at site 
2. For fast monomer diffusion, we estimated the reorganization time in the main text as half the average material 
age. 

When diffusion is slow, we must estimate the reorganization time differently. For equilibrium polymers, after the 
switch in nucleation, site 1 comes to a local quasi-steady-state free monomer concentration: as diffusion slowly drains 
the free monomer pool at site 1, depolymerization keeps the free monomer concentration roughly constant. The 
quasi-steady-state free monomer concentration Ci is 

ci = ^ . (SIO) 

k+ 

At site 2, the free monomer concentration also quickly comes to a value that stays roughly constant for a long time. 
However, this value is smaller than ci since nucleation removes free monomers from the pool at site 2. We may 
estimate this free monomer concentration C2 roughly by 

k-N2 - k+C2N2 - 2knucA = 0, (Sll) 
fcnucici-^^^ = 0, (S12) 

where A^2 is the number of polymers at site 2. Although A*2, like C2, is time-dependent and slowly grows as diffusion 
moves monomers from site 1 to site 2, immediately after the switch A^2 quickly comes to a value that stays roughly 
constant for a long time. Using this value, which is approximately half the final expected steady-state number of 
polymers, and using the final expected MFPT (which is the same as before the switch in nucleation), we may come 
up with a rough estimate for C2 . Then an estimate for the reorganization time is just the time needed for half the 
final steady-state amount of polymerized material to diffuse down the concentration "gradient" from site 1 to site 2. 
If 

Csteady is the final steady-state concentration of poljrmerized material, then our estimate for the reorganization time 
may be written as 



^steady/ 2 
~ kD{ci - C2) 



cquil 



(S13) 



For typical parameter values (such as those used in Fig. 2D), we find that this estimate for the equilibrium-polymer 
reorganization time is accurate to within a factor of 2. 

For active polymers, after the switch in nucleation, site 1 can replenish the free monomer pool only by polymer 
disassembly. Since diffusion is slow, monomers released by disassembly typically rejoin other polymers at site 1. Hence 
after the switch, the number of active polymers at site 1 rapidly drops while the average polymer length grows, until 
eventually site 1 has only a few very long, very long-lived polymers. These few long-lived polymers then quickly 
accumulate and hydrolyze all the free monomers at site 1, and then all the polymers disassemble. Therefore, a rough 
estimate for the reorganization time is just the time needed to hydrolyze all the monomers at site 1, plus the time 
needed for half of that material to diffuse to site 2: 



^active ^ ^total,! . 1 tQT A\ 

where Ctotai.i is the initial total concentration of material at site 1. For typical parameter values (such as those used 
in Fig. 2D), this estimate for the active-polymer reorganization time is accurate to better than 10%. 



COMPARISON OF EQUILIBRIUM AND ACTIVE-POLYMER LENGTH DISTRIBUTIONS WITH THE 

SAME AVERAGE LENGTH 

In Fig. 2B, we showed that equilibrium and active polymers with the same MFPT have qualitatively different length 
distributions: the active polymers have a much longer average length. Figure S2 shows that even for equilibrium and 
active polymers with the same average length, there are more very long active polymers than very long equilibrium 
polymers. 
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Figure S2: Histogram of lengths for equilibrium (blue) and active (red) polymers with average length ~ 90. Black curves are 
theoretical fits from (5) and (7). 
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Figure S3: MFPT versus average length for equilibrium (blue) and active (red) polymers, with the steady-state fraction of 
polymerized material held at 95%. Curves are theoretical fits using MFPT — {L)^^^^^ for equilibrium polymers and (6) for 
active polymers. 

MFPT VERSUS AVERAGE LENGTH FOR 95% POLYMERIZED MATERIAL 

In the simulations described in the main text, we adjusted the nucleation rate fcnuci so that the steady-state fraction 
of polymerized material was always held fixed at 75%. For comparison. Fig. S3 shows the MFPT of equilibrium and 
active polymers as functions of their average length, akin to Fig. 2A, with the steady-state fraction of polymerized 
material held instead at 95%. The general trends are the same as those seen in Fig. 2A: for a given average length, 
active polymers disassemble significantly faster than equilibrium polymers in a manner well-predicted by Fokker- 
Planck theory, and the difference increases as average length increases. The MFPT is measured in the same time 
units as defined in the main text, even though the steady-state free monomer concentration is different. 

MFPT VERSUS AVERAGE LENGTH FOR PARTIALLY- VECTORIAL HYDROLYSIS AND FINITE 

DEPOLYMERIZATION VELOCITY 

In the simulations described in the main text, we studied a minimal microtubule model with vectorial hydrolysis (i.e. 
only the last monomer in the polymer cap is subject to hydrolysis) and infinite depolymerization velocity when the 
polymer cap size goes to zero. Here we relax these assumptions. First, we consider the more realistic partially-vectorial 
microtubule model described by Flyvbjerg, Holy, and Leibler (FHL) [1]. This model has a vectorial hydrolysis rate kh 
for any GTP-bound monomer adjacent to a GDP-bound monomer (i.e. identical to our minimal model), and a smaller 
non- vectorial hydrolysis rate r = kr/Sx for all other GTP-bound monomers, representing the rate per unit length at 
which new GTP-GDP interfaces are created. 6x ~ 8nm is the length scale of an individual tubulin monomer. In all 
other respects the FHL model is identical to our own. 
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Figure S4: MFPT versus average length for the FHL model [1]. Curve is a numerical fit to a power law, MFPT ~ (L) 
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Figure S5: MFPT versus average length for the modified FHL model [1] with finite zero-cap depolymerization rate. Curve is a 
numerical fit to a power law, MFPT ~ ■ 

Figure S4 shows the MFPT versus average length for polymers in this more reahstic microtubule model, akin to 
Fig. 2A. In these simulations kr — fc/j/lO, roughly the value FHL found by fitting to experimental data. The MFPT 
scales ^ (L)*^'^, in between the ~ (L)^^^ scaling of our purely vectorial model and the ^ {L) equilibrium scaling. 

Since the FHL paper did not consider microtubule disassembly dynamics (the purpose of the paper was only to 
predict the distribution of polymer catastrophe times), we next consider the effect of adding to this model a finite 
depolymerization rate fc!^'''''* representing the rate of polymer disassembly after the cap size goes to zero. We choose 
^(dep) twice the polymerization rate k^c. All other parameters are the same as in the other FHL simulations 

above. Figure S5 shows the MFPT versus average length in this modified FHL model. Again the MFPT scales 

Hence, the MFPT scales sublinearly with mean polymer length in a more realistic non-equilibrium microtubule 
model. We did find two regimes that produce linear scaling in the modified FHL model: if the average depolymerization 
time for the entire polymer is much longer than the average time for the cap size to reach zero, then the MFPT is 
dominated by the depolymerization time and scales linearly with average length (data not shown). However, this is 
not the regime observed experimentally for microtubules [2]. Also, if the vectorial hydrolysis rate goes to zero (i.e. 
the hydrolysis rate is the same for all GTP-bound monomers) , then the MFPT scales linearly with the average length 
(data not shown), but this is not thought to be the case for microtubules. 
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Figure S6: MFPT versus average length for the two-species equiUbrium model in the regime where one monomer type is 
energetically favored over another. This regime is designed to mimic our non-equilibrium model from the main text. Curve is 
a numerical fit to a linear relationship, MFPT ~ (L) . 



MFPT VERSUS AVERAGE LENGTH FOR TWO-SPECIES EQUILIBRIUM POLYMER 

In the main text we argued that a large class of equilibrium polymer models produces mean lifetimes that scale 
linearly with mean lengths. This class includes all equilibrium polymers comprised of an ordered sequence of monomers, 
each of one of q types or conformations (or q'^ types for k protofilaments), where interactions between neighboring 
monomers {i, i + 1} contribute free energy J^i ^^i^ to the total free energy of the polymer. This class maps one-to-one 
onto 1-d, g'^-state Potts models, and produces a free energy that scales linearly with mean length. 

To provide further validation for our argument, we simulated one such equilibrium model with q = 2 different 
monomer types. This model is designed to be the equivalent of our minimal non-equilibrium microtubule model, but 
subject to the equilibrium constraint of reversibility. Specifically, monomers of one type (the "GTP-bound" type) 
can attach and detach at the front of the polymer. Monomers at the back of the cap can interconvert between the 
two types, and when the cap size goes to zero (i.e. when the polymer is entirely composed of monomers of the 
"GDP-bound" type), monomers of the "GDP-bound" type can both attach and detach. 

First we consider an energetic regime comparable to our non-equilibrium model. We set the energy of "GTP- 
bound" monomers to favor attachment, and that of "GDP-bound" monomers to favor detachment. This means 
that "GTP-bound" monomers are more stable in polymers, and "GDP-bound" monomers are more stable as free 
monomers. However, in order for monomer interconversion to occur mainly from "GTP-bound" to "GDP-bound", 
as for the non-equilibrium system, "GDP-bound" monomers must also be more stable than "GTP-bound" monomers 
even in polymers. This implies that the total monomer population will consist mostly of "GDP-bound" monomers. 
Indeed, when simulated in this regime, with the same overall and fraction of polymerized material as in our 
original simulations, we find a ratio of "GDP-bound" monomers to "GTP-bound" monomers of 10:1. Furthermore, in 
the limit where monomer interconversion approaches unidirectional from "GTP-bound" to "GDP-bouird", the system 
approaches a single species, and we recover our original equilibrium model. Finally, Fig. S6 shows the MFPT versus 
average length for this two-species equilibrium model. The MFPT scales linearly with average length as expected. 

Second we consider the regime where the two monomer species are energetically equivalent. Figure S7 shows the 
MFPT versus average length for this model. Again the MFPT scales linearly with average length as expected. 
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Figure S7: MFPT versus average length for the two-species equilibrium model in the regime where the monomer types are 
energetically equivalent. Curve is a numerical fit to a linear relationship, MFPT ~ {L) . 



